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1. Introduction 

Our knowledge of a, is of fundamental importance. This statement, which is obviously true 
in QCD, can potentially go far beyond the Standard Model. Let us remember, for instance, that 
the value of as{Mz) strongly affects the physics at much higher scales in the context of Grand 
Unification[|l|]. An increasingly accurate determination of is, therefore, an important goal which 
fully deserves the continuous effort it has attracted [^]. 

When determining a^, the decay of the tau lepton has proved to be a very convenient process 
[^]. For one thing, it is the only lepton which is heavy enough to decay into hadrons, so its decay 
obviously happens without the complications of a hadronic contamination in the initial state. Fur- 
thermore, thanks to a tremendous experimental effort, there are also very accurate data available 

The flip side of all the previous advantages is that the tau mass, ~ 1 .77 GeV, is very close to 
the QCD scale, Aqcd ~ 1 GeV, which means that effects going beyond perturbation theory have to 
be brought under very good theoretical control in order to be able to reliably obtain a very accurate 
determination of a,. Or else, discrepancies arise. For instance, recent determinations of find the 
following values: 

ttsiml) = 0.344 ± 0.005exp ± 0.007th [5] , a, (m^) = 0.332 ± 0.005exp ± 0.015th [6] , 
a,(m^) = 0.321 ±0.005exp± 0.012th [7] , a,(m2) = 0.316 ±0.003exp ± 0.005th [8] . (1.1) 



With the high level of precision claimed, these determinations are not compatible. So, either there 
are assumptions which went into the corresponding analysis which are not correct, or the errors 
have been underestimated. 

Perturbative effects (even after resummations) are not all there is. After all hadrons cannot 
be obtained solely from summing Feynman diagrams and there are also genuinely nonperturbative 
effects. These are understood mostly in the framework of the Operator Product Expansion (OPE), 
where they take the form of condensates [^[10], but there may be others[pl|]. In this framework 
there is a potential problem since the OPE is supposed to be a valid expansion only in the euclidean 
(i.e. > 0), whereas the experimental data are fully restricted to the minkwoski region (i.e. 

< 0). In order to connect the two regions an assumption about analytic continuation has to be 
made. For inclusive quantities like tau decay, it amounts to using Cauchy's theorem on the OPE 
instead of using it on the full Green's function, which is the one truly satisfying the right analytic 



properties. This assumption, first proposed in Ref. [12], has been employed in all analyses of 
tau decay up to date, starting with the pioneering work in Ref. The amount by which this 
assumption fails is referred to as Duality Violations (DVs) (see Eq. ( |l.4| ) and below for a more 
precise definition). Obviously, if this assumption is not right with enough accuracy, an associated 
theoretical error should be included in the final error for a,. 

It is clear that the properties of the OPE are crucial for understanding the physics of tau decay. 
For instance, a recent analysis based on Aleph data finds for the gluon condensate the following 
values: 



^{GG) 



= (-0.8 ± 0.4) X GeV^ , — (GG) 

Vector TT 



= (-2.2 ±0.4) X 10^2 GeV^ , (1.2) 

Axial 
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depending on whether the vector or the axial-vector spectral functions are employed. These two 
values are not compatible and, if taken literally, they would signal a clear breakdown of the OPE. 

Given the situation, and the precision sought, we think it is time for a reassessment of the total 
theoretical error involved. Studies have been made of the error involved in the different resumma- 
tions of perturbation theory [^, |^], and also of the combination of spectral functions which may be 
optimal [^, [ij, but the amount of work dedicated to assess the validity of the assumption that DVs 
are absent is surprisingly meager. This is so even though there are rather general expectations that 



it must fail at some level [15]. This is partly why we decided to take a fresh look at this problem 



and make a reasonable estimate of the theoretical error involved 1 16 - 18]. 
Defining 

the central equation for the discussion is given by [|l^, |T8|] 



'0 . .si „ . . 1 



dsP{s) -imnv,^(^) = -— _^ j^^p(^^)n^^^(^^) + ^|r;(^o) • (1-4) 



=«0 



where 



(5o) = - / dsP{s) -ImAvAs) . (1.5) 

is the function encoding all the DVs. In particular, if the OPE were a convergent expansion in 
the region of interest, one would have AvAl^) = and no duality violations. Effectively, this is 
the assumption made in all analyses up to now. The problem with this assumption is that, on the 
minkowski axis, where the data are located, this obviously cannot happen as the OPE does not 
reproduce the spectrum, i.e. we know that ImAvAil^) 7^ 0. Therefore, the question to address is 
rather how large ^^^(^o) can be. Regretfully, and this probably explains why the question has not 
been addressed before, there is no theory of DVs so the answer cannot be obtained in QCD from 
first principles. 



2. Extracting Duality Violations from z data 



The function ImAy^(i^) in Eq. (L5) requires knowledge of the spectrum up to infinite energy 
and this is, in principle, unknown. This is the main difficulty. In order to make progress, it is 
necessary to have a physically motivated model with which to extrapolate the data from the region 
below the tau mass to the region above it. 

Motivated by very general arguments on Regge theory and asymptoticity of the OPE, we would 
like to suggest the following simple form for the parametrization of the DVs: 

-ImAv,A(^) = e{s- s,ni„) e^^^' sin {av,A + I^va^) ■ (2- 1) 

71 

The exponential fall-off is expected to originate from the intrinsic error inherent to an asymptotic 
expansion and the sine function from a harmonic expansion of the periodic function representing 
the daughter repetitions in the spectrum of Regge theory. The step function ensures the validity of 



the parametrization (p7l|) only at high-enough energies. In other words, (2.1) represents the first 



correction to the asymptotic behavior of the spectral function at high energies given by the parton 
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Figure 1: Different combinations of spectral functions (V, A, and V ± A) (data points in black) compared 
with Fixed-Order perturbation theory (green flat hne) and the result of our fit (2. 1 -2.2) (red curve). 



model plus condensates. In Refs. [15|-|18] a model is studied which realizes all the expected 



features known of the OPE in QCD and Regge theory, and for which Eq. (gjj) is the correct 
behavior. 

A fit of perturbation theory^ plus the parametrization (2J) to the spectral functions, in the 
window 1.1 GeV^ < ^ < m^, yields [16|: 



K-y = 0.018 ±0.004 
TV =0.15 ±0.15 GeV" 
ay =2.2 ±0.3 
j8v = 2.0±0.1 GeV-2 
10 



dof 79 



0.13 



Ka = 0.20 ±0.06 

Ya = 1.7 ±0.2 GeV-2 , 

a^ = -0.4±0.1 

j8A = -3.0±0.1 GeV-2 



dof 



17 
78 



:0.22 



(2.2) 



The result of this fit is compared to the data in Fig. |TJ As one can see from the plots and the 
corresponding x^/dof, the quality of the fits is more than acceptable. 

Since now one has a way to extrapolate the spectral function all the way to infinity, it makes 



sense to ask whether the ansatz (|2JJ) together with the values for the parameters ( |2.2[ ) satisfies 
the 1st Weinberg sum rule^ (WSRl).^ One may not expect, however, that the condition to satisfy 
WSRl will do away with DVs altogether. The reason is that a, is determined from the V ± A 
combination of spectral functions in T decay whereas WSRl obviously depends on the orthogonal 
combination V — A which contains no perturbative contribution (and, therefore, no dependence on 



'Condensates give negligible contributions [jl^. 

^We remark that the 2nd Weinberg sum rule breaks down away from the chiral limit. 



We thank A. Pich for asking this question. 
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Figure 2: Left panel: result of the WSRl for sq — mj. employing T data. Right panel: the WSRl extrapolated 
to higher energies with the use of our ansatz (2. 1-2.2). The (blue) horizontal line is given by the value of 



Us). At any rate, due to the smallness of the u,d quark masses, the WSRl may be expressed, within 
a very good approximation, as 

1 

271 



WSRl(5o) 



ds Imnv_A = 



(2.3) 



This sum rule, which is supposed to be valid only for very large, follows from the general result 
in Eq. (1.4), provided the DV contribution ^fz^i^o) is set to zero.^ This can be seen by taking 
the particular polynomial P{q^) = 1 and remembering that there is no condensate of dimension two 
from the OPE in the V — A combination. In other words , for any finite sq, the WSRl is a measure 
of DVs at that scale since only when 



oo can one make sure that DVs vanish (see Eq. (1.5)). 



The result of the integral up to an scale is shown in Fig. § On the left panel, one sees that the 
experimental data grossly violates the WSRl at s = m^, exposing thereby the existence of DVs at 
this scale. At higher energies, the assumption of no DVs is tantamount to essentially a constant 
zero line from onwards for the V — A combination of spectral functions shown on the rightmost 
lower panel in Fig. This means that the result obtained on the left panel of Fig. ^ for = 
cannot improve at higher energies but remains constant. Our conclusion, therefore, is that DVs are 
clearly not vanishing at m^. 

On the contrary, when the data are extrapolated at higher energies taking into account DVs 

), the sum rule does get satisfied within errors, as it should. This 



with our parametrization (|2. ll^^ 
is shown on the right panel of Fig. |2[ 

The presence of the DV term ( |2.l| ) with the values for the parameters extracted from the fit 



( |2.2| ) affects the standard extraction of made through Eq. ( |1.4| ) because of the contribution 



coming from ^v^i^o)- We have estimated in Ref. Qlq ] the associated theoretical error in due to 
DVs as 



da, 



nir 



0.003 - 0.010 



(2.4) 



where the spread of values includes the propagation of all the errors involved. Another attempt 
at estimating this theoretical error was made in Ref. [^], with the result that DVs were negligibly 
small. However, unlike in our analysis, no detailed fits to the spectral data were made in this 
reference. 



Contributions from condensates are expected to be negligible. 
'Again, the contributions from the condensates are negligible nil 
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Figure 3: Plots of the vector spectral function in tau decay (left panel) and e+e (right panel) as compared 
to perturbation theory (green flat line) and the result of the combined fit (3.1) (red oscillating curve). 

3. Inclusion of e+e data. Conclusions 

Since our DV ansatz allows one to go to higher energies, one may wonder how well our vector 
spectral function compares with the data extracted from e^e^.^ Of course, our ansatz ( ^!1| ) cannot 
be directly used for e^e^ data because it contains only isospin one (i.e. the ud combination) 
whereas in e^e^ one has the flavor singlet combination {2/3)uu — {l/3)dd — (l/3)ss. Neglecting 
this fact and adjusting only for the different charges involved (a bold step to take) plus a shift in the 
parameter a — > a' (to take into account that the s quark is much heavier than the u, d; a modification 



which is suggested by the model in [15]-[18]) we have performed a simultaneous fit to the vector 
spectral function from T decay, together with the vector spectral function in e+e^ above 4 GeV^ 
up to the charm threshold [^. The e^e^ data between and 4 GeV^ is controversial [21] and 
we have not used it. 

The result of this simultaneous fit becomes 



Jv 
ay 

dof 



0.024 ±0.004 

0.40 ±0.12 

1.82±0.19 

2.14±0.11 

5.2±1.4 

22 

— -0.25 , 
87 



(3.1) 



which entails a shift from the values obtained only with T data in ( |2.2| ), although it is compatible 
within errors. The result of the fit (3J) can be seen in Fig. ||. As one can see, given the assumptions 



made, our ansatz is not grossly incompatible with e^e data. Taking the values (3J) at face value, 
the theoretical error in turns out to be smaller than (2.4) by a factor of 3. 



We think it is advisable not to fall on the optimistic side when it comes to estimating a theoret- 
ical error. From the plots in Fig. |T] and ^ we conclude that a theoretical error in the determination 



of as from duality violations at the level of Eq. (2.4) is not at all excluded. However, this error has 
not been included up to now in any determination of a, from T decay. Until we learn more about 



DVs, we think it should. We refer to Ref. [ 16] for more details and further discussions 



Regretfully, there are no data at higher energies in the axial channel. 
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